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CHAPTER I 
INTRODUCTION 

In 1916, Albert Einstein published his famed thoery of general rela- 
tivity. Although many contemporary physicists and astronomers based their 
work on the principles of this theory, proof of Einstein’s description of 
the universe remains inconclusive. In fact, alternatives to Einstein's 
theory, most notably the Brans-Dicke scalar-tensor theory, have been deve- 
loped.^ 

The Gravity Probe experiment, conceived by the late Professor Leonard 
Schiff of Stanford University, seeks to test Einstein's theory, and has 
been recognized by the U. S. Government as having far-reaching scientific 
importance. Since 1962, the government has supported the development of 
this experiment in which academic, industrial, and governmental institu- 
tions have participated. Currently, NASA's Marshall Space Flight Center is 
charged with the overall responsibil ity of the experiment. 

In proposing this experiment. Professor Schiff had shown that a gyro- 
scope in orbit undergoes a relativistic precession in the framework of the 
2 

fixed stars. In the "Gravity Probe B" (GP-B) satellite, four almost 
perfectly spherical fused-quartz gyros, each the size of a ping-pong ball, 
will be coated with niobium, and element made superconductive at tem- 
peratures near absolute zero. In this condition, the gyros can be suspended 
electrostatically in the weightlessness of space. It is expected that the 
drift rate of these gyros will be held to less than .001 arc-second a year. 

For gyroscopes in polar orbit, there are two predicted precessional effects 
which are perpendicular to one another, as shown in Figure I-l. The larger 
effect is due to the motion of the gyroscope through the Earth's gravitational 
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field and is referred to as geodetic precession. The other predicted 
effect is the spin-spin coupling between the gyroscope and the Earth's ro- 
tation and is referred to as the motional effect. As shown in the figure, 
the geodetic effect is in a plane of orbit approximately parallel to the 
Earth's axis, The motional effect lies in a plane perpendicular to both 
the orbit plane and the Earth's spin axis. The predicted magnitudes of 
procession vary ffom theory to theory. The purpose of this experiment is 
to determine the most valid theory by comparing the experimentally 
measured magnitudes with the various magnitudes obtained from different 
predictions. 

The most important part of the GP-B satellite is a dewar containing 
the four spherical gyroscopes, a telescope which tracks the star Rigel, 
and an amount of supercooled helium. Spacecraft attitude and trans- 
lational control is accomplished by proportional thrusters with the pro- 
pellant supplied by the boil off helium from the dewar. The control 
signals are derived from the telescope when the guide star is visible, and 
frcm the gyroscopes when the star is occulted 

The experiment design goal is to make drift measurements accurate to 
.001 arcsecond per year, an unprecedented requirement for drift measure- 
ment. To ensure the success of the experiment, every conceivable error 
source in the GP-B system needs to be identified and its effect on drift 
measurement investigated. 

In April 1980, NASA/MSFC awarded a study contract, NAS8-33849, to the 
University of Tennessee to initiate a system approach to GP-B error analysis. 
The study results included a block diagram error model, identification and . 
discussion of source errors, some preliminary error analysis, documentation 
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of error budget, and several recommendations. These results provided an 
organized picture of various errors Involved. 

The present contract, NAS3-34426, Is almed at an 1n-depth Investiga- 
tion of Individual source errors and their effect on the accuracy of the 
GP-B experiment. This report documents the results of the study performed 
over a fourteen month period, from 24 June 1981 to 31 August 1982. The 
main emphases during this report period were placed in the following three 
areas; 

1. Refinement on source error identification and classifications 
of source errors according to their physical nature. The 
result is given in Chapter II of this report. 

2. Error analysis for the GP-B data processing. The result 
obtained 1s given in Chapter III. 

3. Measurement geometry for the GP-B experiment, to be in 
Chapter VI. 

It should be mentioned that this 1s a progress report, therefore the 
results reported can be further Improved. 

References 

1. C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation , W. H. 
Freeman and Co., 1973. 

2. L. I. Schiff, "Mdtlon of a Gyroscope According to Einstein’s 
Theory of Gravitation," Proc. N.A.S. , VoT. 46, 1960, 871-882. 

3. 0, C. Hung, "Gravity Probe B Error Analysis," Final Report 
for NASA Contract NAS8-33849, The University of Tennessee, 
Knoxville 18 February 1981. 
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CHAPTER II 

PROPERTY AND CLASSIFICATION OF 
GP-ERROR SOURCES 

] 

Figure II-l Is a block diagram error model for the GP-B experiment. 

The model was developed in the previous contract based on the data con- 
tained In Reference 2, This diagram gives a "macro" description of errors 
where error sources are lumped into nine groups entering the system at 
nine entry point--. During this contract report period, further effort 
was exerted on error source identification using the data available in 
Reference 3. In addition! these error sources were classified in two 
different ways based on their physical properties. The purpose of 
classification Is to expose the possibility of suppressing the effects of 
individual errors and to suggest techniques for such suppression. 

The errors can be classified according to there statistical nature 
as follows; 

1, Deterministic type — This type of errors can be compensated. 

2, Random variable type — This type of errors are uncertain 
constants. The effect of some, but not all, of these errors 
can be eliminated by rolling the spacecraft and by orbit 
averaging. 

3, Random process type — This type of errors aj'e uncertain 
time-varying quantities. Some of their effect can be re- 
duced by averaging or by filtering. 

The errors can also be classified according to their physical forms in 
the following way: 

1. Bias — The effect of some of this type of errors may be 
reduced by roll averaging and orbit averaging. 
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Figure II-l. A block diagram error irodel. 








2. DHft. 

3. Scale factor uncertainty — May be Identified by Kalman 
filtering, 

4. Noise — May be helped by Kalman filtering. 

Table II-1 Is a list of error sources, their classification, and 
methods for the reduction of error effect. In the table, r.v. means 
random variable, r.p. means random process, 1/f means 1/f noise, and det 
means determlnistics. Notice the estimated errors after data processing, 
their root-sum-square (RSS) total Is 1.253 mllliarcseconds assuming zero 
proper motion for Rigel. It is expected that future progress In hardware 
and software development will reduce the RSS error. Notice also that 
the table Is open-ended, vith blanks, “TBD", and question marks to be 
taken care of by further study. 

Figure II-2 shows the budget curve for the average individu.i error 
as a function of the number of equal size error sources. The curve is 
based on a RSS error of 1 millarcsecond. For example, if the total number 
of equal size error source is 44, the individual error source should not 
cause more than .15 milliarcsecond error for the GP-B experiment. 

References 

1. 0. C. Hung, "Gravity Probe B Error Analysis," Final Report for 
NASA Contract NAS8-33849, University of Tennessee, Knoxville, 

18 February 1981 . 

2. 0. R, Parker, "Gravity Probe B Phase A Report," National 
Aeronautics and Space Administration, Marshall Space Flight 
Center, March 1980. 

3. C. W. F. Everitt, "Report on a Program to Develop a Gyro Test 
of General Relativity in a Satellite and Associated Control 
Technology," Stanford University, Hansen Laboratory, June 1980. 
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6. D, H. Sheingold, Analqq-digital Conversion Handbook , Analog 

Devices, Inc., Norwood, Mass, 197^. ^ 

7, L. Motz and A. Duveen, Essentials of Astronomy , Columbia Uni-* 

versity Press, New York, 1971. ' 
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CHAPTER in 

FINITE-WORDLENGTH INDUCED ERRORS IN 
KALMAN FILTERING COMPUTATION 


1. Introduction 

The problem of finite-wordlength effect on digital computations has 
been investigated extensively during the past twenty years. Finite- 
wordlength property of a computer requires either rounding or chopping to 
be used to limit the wordlength of a number. Since most computers use 
rounding technique, only rounding errors will be considered in the sequel. 

There are two approaches to analyze the rounding error, the first 
approach considers the statistical nature of rounding errors, and treats 
them as noise generated in the system, This approach has been widely 
used by those in the field of digital signal processing. In the statis- 
tical error analysis one is usually after the ensemble average and stan- 
dard deviation of the final error based on the estimated characteristics 
of source errors and their propagation through computation steps. This 
approach does not seem to be sufficiently reliable for the analysis of 
GP-B data reduction errors for two reasons. First, GP-B's four experi- 
ment gyros represent only a small sample, their combined statistical 
characteristics may deviate a good deal from those of the population 
statistics, Thus the use of sStatistical analysis here may not give a 
reliable result, Secondly, the GP-B data reduction involves Kalman 
filtering and other rather complex computations. The exact statistical 
nature of rounding error generation by and propagation through these 
computations is not easy to establish. Therefore a more conservative 
approach Is needed. 
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The second approach 1s to establish bounds for the rounding errors 
Involved In computation. This approach provides a very conservative, 
though rather pessimistic, result for rounding error analysis. This 
approach has often been used by those doing numerical analysis, Be- 
cause of the unusual precision required of the 6P-B and the expensiveness 
of the experiment the use of error bound approach provide a much more 
reliable results for the error analysis. Therefore this approach will 
be used for ensuing rounding error analysis. Since Kalman filtering 1s 
the main activity In GP-B data reduction, the present chapter 1s devoted 
to the analysis of rounding error In Kalman filtering computation. 


2. Rounding Procedure In Floating Point Representation 
Let X be a number 

X = (t'd-jdg — ) X b® (1) 

v/here b is the base of the number system used and e, an integer, Is the 
exponent. In general the mantissa part of the number may have infinite 
number of digits for an exact representation, such as for The number 
(1) may also be represented in the form 

X = u-b® + vb®"'^ (2) 

where ^ |u| < 1, 0 < |u| < 1, and u contains only t digits. 

Examples; Base 10 numbers. 

(a) 12.3456 = .1234 x 10^ + .56 x 10'^ 

Here b = 10, ,t = 4, and 3 = 2 

(b) -.0123456 = -.1234 x 10'^ + (-.56) x 10"® 

Here b = 10, t = 4, and 3 = -1 

The rounding procedure drops off the second term on the right side 
of (2) by appropriately adjusting the value of the first term. 'Thus, 
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after rounding, x becomes X|^ which has a t-digit mantissa *did 2 * d^ 
and an exponent b®. The conventional round-off procedure for any number 
Is as follows: 


r u.b^ 


X = 


u b® + 


I U b® - b®-^ 


1f |v| < i- 
if V >i 

1f Vl-1 


(3) 


Note that u and v always have the same sign. 

Examples: b = 10 and t = 4 

(a) X = 765.4567 = .7654 x 10^ + .567 x 10"'' 

Here v ^ J- and e = 3, so 

X = u b® + b®"'^ = .7654 X 10^ + 10^"'' = .7655 x 10^ 

(b) X = 123.426 = .1234 x 10^ + .26 X 10"' 

1 

Here v < p so 

X = u b® = .1234 X 10"^ 

(c) X = - 765.4567 = -.765 x 10^ -.567 x 10"' 

I 

Here v _< - ■g'* so 

X = u b® - b®'* = -.7654 X 10^ - 10^"'' = -.7655 x 10^ 

These results are intuitively obvious. The reason for going through the 
formulations of Equations (1), (2) and (3) Is to pave a way for the sub- 
sequent analysis of rounding errors. 
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3. Rounding Errors 1n Floatinc? Pointing Representation 
The "absolute rounding error" In Is defined as 

|x|= > 0 (4) 

From (2) and (3), It Is clear that 

(5) 

Examining (1) shows that |u b®| ^ b®"^ because u ^b"^; and |x| ^ |u b®| 
because the second term, having similar sign, Is dropped. Hence 


DEfine the "absolute relative rounding error" e as 


G = 


Xr-x 


_ Ixl 

-'nr 


By (5), (6) gives 


G 5, ^ b"^ ^ ” B 


( 7 ) 

(8) 


The quantity p 1s called the "unit rounding error" which represents the 
absolute bound of rounding error in the floating point representation of a 
number of base b and having a t-digit mantissa. It is an important 
parameter in the analysis of rounding errors. 

Example: Consider b 10 and t = 4 

Then B 

Let X = 767.4567 = .7654 x 10^ + .567 x 10"’ 

then X = .7655 x 10^ 


IX) = |Xj^-x| == .0433 
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760i7 = -56568 X 10-^ < e 

For the sake of comparison, the chopping error In floating point 
representation of a number will be analyzed next. 


4. Chopping Error 1n Floating Point 'Representation 

For a floating point number In the form of (2), a t-digit chopped 
number Is given by 

= u b® (9) 

Define the "absolute chopping error" x as 

w 

ix:j= ivxi = M®'* 00) 

Since |x| <1 

01 ) 


Define the "absolute relative chopping error" as 

l^cl 

-R- 

c1 early, 


( 12 ) 



where is called the "unit chopping error, 
shows 


03) 

Comparing (13) and (8) 


= 26 , (14) 

Example: b = 10 and 6 = 4 
Then 6;, = 10^'^ = 10"^ 

Let X = 765.4567 = .7654 x 10^ + .567 x 10'^ 
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then Xp « .7654 x 10“^ 

V 

IXjl» Ix^-xl = .567 X 10"'' 


5. Rounding Error in Basic Computer Arlthmatic Operations 

For the convenience of the subsequent analysis, notation for rounded 
floating foint number is defined here in two equivalent forms. Let x be a 
floating point number. The rounded value of x is denoted by x or f)i(x). 

Let denote any of the four basic arithmatic operations +, x, and 
/ . The computer value of x*y is f£(x*y), which is related to the exact 
value x*y by 

fMK.^y) = (x*y)(l+e) (15) 

where e is the actual relative rounding error. The absolute relative error 
in (x*y) is bounded by 


e - 


fg>(x*y) - (x^y) 
(x*y) 


< 0 


where E is the unit rounding error. 


(16) 


6 . Rounding Error in Composite Computer Arithmatic Operations 
Repeated Additions and subtractions consider the sum 
s = x^ + Xg + X 3 + x^ 

= ((x^+Xg) + X3) + X4 

The rounded value is 

s " {[(x^+X2)(l+e-j ) + X3] (l+Eg) ^4^ 

= (x-j‘t-X2)(l+e-j )(l+E2)(l+£3) + X3(l+e2)(l+e2) + X4(l+&3) 
ck (Xi+X2)(l+ei+e2tE3) + X3(l+£2+e3) + x^(l+£3) 

18 
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The rounding error is 

S = (XiX2)(e^+e2*^e3) + X3 (62+23) t (63) 

« (X1+X2+X3+X4) (6^+62+62) - X3E^ - X^ (e-|+e2) 

The absolute relative rounding error is founded by 
s 1 33 " 3X3 *• 20x4 ^ 33 

In general , for a sum of n terms 
n 

s = I X . 
j=l ^ 

the absolute relative error is bounded by 
e ± (n-1 )3 


(17) 


(18) 


Repeated Multiplication and Division . Consider the following 
combination of product and quotient 




^ 1^2 
^1 “ 




the rounded value is 


Q == 




(i+n-|)2i“"^ (1+E^+Ti^) 


where and are relative rounding errors due to multiplication and 
division, respectively. The rounding error and the absolute relative 
rounding error are, respectively, 
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For the general case 
X^Xg — - X 


Q = 


n 


y^yg — y^ 

The absolute relative sounding error is bounded by 
p < (n+m-1 )3 


09 ) 


( 20 ) 


7» Norms of Vectors and Matrices 

Norms of vectors and matrices are useful In the analysis of rounding 
errors in matrix operations. • The following definitions of norm will be 
adopted in this study, 

For an n-vector ^with elements Xj» define the vector norm as 
I |x| I “ Max|xJ (21) 


Clearly, the norm has the following properties: 

(1) ll>ill>0 

(^1) Hill 1 “0 oiiiy if X “ 0 
(111) ||x+y|| < ||x|| + I [ill I 
(1v) I lax.! I ® l®l’l|x|l for any a eR 

For a maxn matrix A with elements a^.^, define the matrix norm as 


n 

I 

j i^l 

This norm has the following properties: 


|A|| = flax I Ja^j 

i -f-l ‘J 


( 22 ) 


(I) l|A|| iO 

(II) 1|A|| =0 only If A = 0 

(III) i|A+B|| < l|Al| + ||B|| 

(iv) I |a A[ I - |a| ■ I lA| |for any aeR 
(V) IIABII < ||A!M|8|| 
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8* Rounding Enron In Matrix Addition 

Let A and B be two nxm matrices, the rounded sum of them Is 
fjt[A+B] « A+B+R (23) 

where R Is the rounding error matrix. By the definition of matrix norm 
(22) and in view of (15), the norm of the rounding error matrix 1s bounded 


IIRII ^3 ||A+B|| (24) 

where 3 “ as given by (8). The relative norm of |1R|| Is bounded 

by 


e 


l|R|| <8 

IIA+Bll - 


(35) 


which Is the same as the relative rounding error of the sum of two numbers 
as shown in (16). 


9. Rounding Error 1n Matrix Multiplication 

Since elements of a matrix product are Inner products of vector 
pairs, the rounding error associated with an inner product will be analyzed 
first. The result will then be used to analyze the rounding error in a 
matrix product. 


9,1. Rounding Error 1n Inner Product 

Consider the Inner product of two 3~vectors a^ and b 

I = - a-j b^ + a^bg + a^b^ 

f 

The rounded value of ^ ]b Is 

i = nc/y 

“ {Ca-|bi(l+e^ ) a2b2(l+E2)](l’l-e3) + a3b3(l+e^)}(l+E5) 
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^a-jb-j (l+ei+e^+eg) + Rgb2 (1+e^'j+eg+eg) + a^bg (1+e^+eg) 

where e^'s are relative rounding errors associated with basic arlthmatic 

A 

operations, The rounding error 1n I Is 

i » a^b-j (e-j+eg+eg) + a^bg (Heg+e^+eg) + a^bg (e^+eg) 

The absolute value of this rounding error Is bounded by 
li| <30 la^b^l + 30 lagbgl f 23 [aghgl 
In general, the aboslute rounding error of the Inner product of two n- 
vectors Is bounded by 

|I| isn{la,b^| + ! (n+2-j) U^bjl} (26) 

d®2 

The expression for the absolute relative rounding error for an 1i 
product appears cumbersome and is not given here. 

9.2, Rounding Error In Matrix Products 

Consider the matrix product C = AB where A is mxn and B Is nxp. 

The number n will be called "interface dimension" for matrices A and B. 
Using the result of (23) the absolute error of thg elements of C is 
bounded by 

lo^jl i 6 (n|a.^|-lb^'^l + nla^gl ’ l'’2j I ("-I )|ai3l 

+ + 2|a^„l.|b„j|) (27) 

Let [c] be a matrix whose elements are [A] be a matrix whose 

elements are|a^j.|, and‘[B] be a matrix whose elements are lb^j|. 
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Then, based on (27), one has 
[c] I 0[A]OCB] 

where the symbol "P‘ means that the comparison Is done on element by 
element basis for the left and Hght hand matrices, and 



Clearly ||d| 1 ~ n* The norm of the rounding error matrix c Is 
bounded by 


4« ^ ^ 

wncr c rur c 


IICII < nellAlHiBlI =^b’-'t||Al|*|!B|| (28) 

Generalize the above result to a product of N matrices 
P ' - - - Mf, (29) 

With Interface dimensions d^, d 2 > '’n-1' 

P, = 

Then the result of (38) Implies the following rounded matrices, with e being the 
worst error, 

^2 = “ MiMgd+d^e) 

and 

P^l = nCM^Mg] = M^M2(l+d,e) 

P 3 = f CPiMj] = P^M 3 n+d 2 c) 

= M,M,M3(l+d^E)(l+d2E) ^ H,M2M3[1+(d,+d2)e] 
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and 


A 

Rounding errors 1n Pj^ Is 

^N ^1^2 - - •* Mj,|(d-j+dg + - - **+dj^^l)e 

The norm of this error matrix Is therefore 


( 30 ) 


(31) 


N«1 N 



The results of (24) and (31) can be used jointly to handle the matrix 
equation containing both products and sums. This will be demonstrated by 
the following two examples. 

Example ! Compute 
R « ABC^D 

where all matrices are nxn. The rounded R 1s 
R « ABC(l+2nB)(ltR) + D(VHc) 

A 

ThtiJ rounding error of R 1s 

R » C(2n+1)ABC + D]e 
and its norm Is bounded by 

||R|| ^0[(2n+l)l|ABC|| + 1|0||] 

Example 2 Compute 
R = ABC + D 

Where A Is nxm, B is mxr, c Is rxs, and D is nxs. Then 

R a ABC[l+(m+r)c](l+G) + D(l+e) 

R « ABC(Hm+r)e + Dg 

||R|| ^BC(l+m+i-)||ABC|| + ||DH] 


24 



ORIGINAL PAGS IS 
OF POOR QUALITY 

These two examles shows that the rounding error norm of matrix 
addition does not Involve the dimension of the matr1ces» but that of 
matrix product Involves all the Interace dimensions, 

10, Rounding Error In Matrix Inversion, First Approach 

-1 

Let A be nonsingular nxn matrix, Its Inverse A satlsffv^s the 
relationship 
-1 

A A “ I, the Identify matrix 

-1 

Let be the jth colume vector of I and h^j be the jth column of A . 

Then hj 1s the solution of 

A^ = Jij j = 1 to n (33) 

Thus A" can be obtained by solving (29) n times using different Uj 
each time. The solution 1s usually done by a method based on the Gaussian 
elimination with partial pivoting. The present concern Is the rounding 
error associated with the conputatlon of A . The analysis will be done In 
two steps: First, find error 1n A computed from the exact A. Second, 

find error In A“^ computed from A’ “ A + AA where aA 1c the error In A. 

' -1 

Rounding Error In A when A Is exact. Let h^j be computer solution 
of (33). Define the "residue" associated with h,- as- 

J 

Hj = fl h. - Uj (34) 

The error In hi| Is 

J 

hj = hj - hj = A-Vj (35) 

The rounding error matrix for the computer inverse of A is 

E = [£l iz - - - h„] (36) 
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Define the “residue matrix" for the computer Inverse 


R ® r« “ 


TIton 


E » A“'’r 

The norm of E Is bounded by 

l|ElUI|A‘''lN|R|| C 

and the relative norm of E 1s bounded by 

8 <I|R(I (- 

-1 

Roundlno Error In A ^ when A-^aA 1$ Inverted. Let A be erred to 
A+Mi then the computer solution \±^ for 


[AHA]x « u, 


j « 1 to n 


must satisfy 


A 

[A+M]^j « jr j d « 1 to n 

where ^ Is the residue, Then 

ij + A’^Mhj ■ PC\^ + A'\rj » hj + A'Vj 
The error In hj Is 

ij • ij “ llj ° ■ fiA hj] (A3) 

Using the notation defined In (36) and (37), (43) gives the error matrix 


E - A'''r - A‘l 


i A-''r - A'^AA'^ 


The norm of the error matrix 1s bounded by 
l|E|| illA-’||-||R|| + l|A-^||^-||/^!l 


the relative error norm 1s 

c * * I l^l I I 1 1 ‘ I UA| 

l|A 'll 


Q 


OfiieiNAL PAGE IS 
OF POOR QUAUTY 

Comparing (44) to (40) shows that the latter is a special case of the 
former where M « 0. Eq.(44) appears elegant, but Its practical useful- 

ness Is 1n doubt. The problem is that the residue matrix R cannot easily 
be obtained. In sidditlon, both (43) and (44) do not explicitly depend on 
any wordlength related parameter, such as, the unit sounding error 3, 


11, Rounding Error in Matrix Inversion, Second Approach 

The usual method of matrix Inversion by a computer Is based on a 
repeated use of Gaussian elimination procedure. The procedure consists of 
two parts, namely, triangularizatlon of a matrix and back substitution. 

The rounding error for each part will be analyzed first, followed by the 
analyze of the resultant error. In the following analysis e will denote 
the wrrst value of any rounding error e^. Th'Xs |e| ^ 

11.1. Rounding Error In Matrix Triangularizatlon 


Consider a 3x3 Matrix Equation 


*11 

*12 *13" 


"xf 


Di 

*21 

*22 *23 


^2 


^2 

*31 

Urn ^ 

*32 *33_ 




>=3 


A matrix 


(45) 


Let a|j(0) = and b^(0) = for i, j = 1 to n. The first step is to 
condition the first column. Let 

( 0 ) 


"'21 

then let 


- 

'll ■ 


“21 

a^jW 


21 


( 0 ) 


aggd) = aggCO) + m2.,a^2(0) ° ° ETf7(^®12^°^ 


( 46 ) 
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Spi (o) 

^22^^^ * fA[a22(0) “* 

9pi (o) 

= *22'°)(Hc) =5^a,2(0) (U3e) 

api (o) 9n>i(0) 

^^22^^^ “ a^”W ■’* ®22^^^ " a^^l6T 

= ag2(t) + a22(l)(3e) - a22(0)(2e) 

Similarly 


323(1) ~ ^23^^^ ^23^-^ “ ^23^^^ 

832(1) = 332(1) + 332(1) ( 3 c) - 832(0) ( 2 c) 

933(1) - 933(1) + 933(1) ( 3 c) - 833(0) ( 2 c) 

A 

After the first step, the error in the new A, designated A(l), is given 

#w 

by A(l) whose elements are 


^ijH) = 


i=l; j=1.2,3 


(47) 


a^j(l) (3e) - a^.j(l) (2e) i,j=2,3 


Define 


A^(l) == 


0 

0 

0 


0 0 
^32^^^ ^33^^^ 


and 


VD = 


0 0 
0 322(1) 

0 332(1) 


0 

823(1 ) 
933(1 ) 
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Then 

A(l) = C^ijO)] = SeA^O) - 2 eAt(0) 
c.3[3At(1) - 2A^(0)] 

Expressing A in terms of A rather than A is important since 
from the computer but not A. 

For the b^- coefficients we have 

a.,(0) 

b^d) = bi(0) = m.ib^(o) = b^(0) -5^b,.(0) 

Then 

b^O) = «[b^n)] = b.(0)(l+e) - 5^|{^bi(0)(l+3e) 

a.i(O) 

= b,-(l) + e b^(0> - 3e 5^|^b,(0) 

= b^(l) + 3e b^(l) - 20 b,(0) 

Define 


“o“ 


~o“ 

b2(0) 

and b-j(O) = 

bgO) 

b3(0) 


bad) 


then 

b(l) = 3e b^O) - 2 g b^(0) 2i 3[3b^(l) - 2 b^(0)] 

A 

Again, expressing ^ in terms of ^ rather than b is important 
available from the computer but ^ is not. 


( 48 ) 

is available 


(49) 

A 

since b is 
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After the first reduction step, one has 


a^l(0) 812(0) 

IS 

CO 


’^1 


b^(0) 

0 

CM 

CM 

< 

0 

8230) 



s 

£2(1) 

1 

0 

0 >> 

CO 

ro 

0 

® 33 ^^^ 

t 



1 cr> 

1 . 


where and are rounded quantities. Jheir errors will be compounded^ 

to the new rounded quantities In the next step of the reduction process. 

The second step of the reduction concerns the second column of the 
matrix. Let 

332 ( 1 ) 

"'32 = - (50) 

322^ I ) 

and let 


K 9 O) . 

333(2) = 833(1 ) + m32 823(1 ) = 833(1 ) - — 823(1 ) 


822 ( 1 ) 


The rounded value Is 


a««(l ) 


833(2) = n[a„(2)] = a„(l)(l+6) - Ji_a„(l)(l+ 3 e) 


' 33 ' 


‘33 


S22(1) ” 


~ ^33(2) + 3 e 933(2) - 2 e 933(1) 


Let 


A 2 O ) = 


00 0 

00 0 

0 0 333 ( 1 ) 
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and 


Ag(2) = 


0 ^33 (2) 


Then 


A(2) = 3bA 2(2) " 2eA2(l) £ii ^[SAgC ) - 2 A 2 O)] 

For the coefficients 1n the second reduction step, 

b3(2) = bgO ) + mgg bjO ) = bgd ) - bgO ) 


The rounded value is 


£3(2) = Mb3(2)] = b3(l)(l+6) - 


^32(1) 


iggO) 

= b3(2) + 3e b3(2 ) - 2b £ 3(1) 


bgdjO+Se) 


Let 



I “1 


0 

bgd ) - 

0 r— 

CO 

and b2(2) = 

0 

b3(2) 


Then 


(Bl) 


1(2) = [3 b2(2)-2 ^d)] B a. g[3 ^(2) - 2 
The resultant errors, 

A = A(l) + A(2) = g[3 A,(l) - 2 A^(0)] + e[3 ^3(2) - 2 Agd)] 
2 


(52) 


= b{ y3A^-(0 - 2A^(i'-l)]}' 


(53) 
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b - io) + b(2) ».3[3b^(l) ^ 2b^(0)3 + 3[3bg{2) . 2^(1)] 




Generalization to an nxn matrix A 


(B4) 


A = A(0) = A(0) = 


®n ■ ■ ■ ®ln 


L^nl ’ " " ^nn-* 


a,j( 0 ) - a^j(O) = 

The matrix obtained after the kth reduction step is 

r 

I I 


A(k) = 




its) 


(56) 


(57) 


A^(k) - 


O CZD 



(58) 


The (n-i)x(n-i) lower right 
diagonal block matrix from A(k) 

The resultant reduction or triangular! zation errors in A and b are, 
respecti vely. 


A = 


C3A,.(i) - 2A^(i-l)]} 


(59) 
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b = sQ C3b^(i) - 2b^(1-l)]} 


(60) 


Finally, the norms of errors due to triangularizatlon are given by 
'n-1 


11*11 - >0' [S|I*|,(H)1| •!|IVW)I|]} 


(61) 


n-1 


l|b|| = [3||b^(k)|| +2|lb^(k-l)||]} (62) 

Note that after the (n-l)th reduction step, the original matrix A has 

A A 

been reduced to an upper triangular form. Denote it by Aj - A(n-l). Thus. 

( 0 ) 

In' ' 



m 


«i.w(k-l) - - a^„(0) 




(63) 


b^ = b(n-l) = [b,(0) bgd) b„(n-r)] 


(64) 


11.2. Rounding Error in Back Substitution 

This problem is approached as follows. Consider the equation 
Ax = b (65) 

where A is an nxn upper triangular matrix. Let ^ be the solution of 

A A 

this equation which contains rounding errors, then find A and ^ 
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such that 

A A A 

A X * b (66) 

A 

has error-free solution x. Thus the rounding error problem has been 
transformed to a error problem caused by perturbations 1n A and b. Define 

AA A-A AX » x-x Ab « t-h (67) 

Then (65) gives 

(A-aA)(x-ax) = £ - a b 

ix = A'^C b- A ^ _ (68) 

In (68), A>L 1s the rounding error in x, A ’ and ^ are computed by the 
computer while solving (65), and A^ and aA are computed from formulas to 
be developed. Assume that A is an upper triangular matrix, the solution 
of (65) involves only the back substitution operations. Rounding error 
due to back substitution will now be analyzed. 


3x3 Triangular Matrix Equation. Consider the equation 

b. 



^12 

a 6 

^13 

0 

^22 

®23 

0 

0 

^33 


I Tvl 
^1 


'1 


(69) 


The equivalent perturbed equation for evaluating rounding errors is 


rr 

a 


11 

0 


*^12 ®13 


®22 ®23 


‘33 J 



A 


— 


<! X 


A 

*’2 


A 


1 

CO 

J 

tmm 




(70) 
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Write (69) as 


^13^3 “ 

^ 22^2 ^ 23^3 “ ^2 

^ 33^3 “ ^3 


From (71c) 


^3 a“ 


The rounded value Is 


’‘3 “ 


' ' "3 ■ 1+F 

Substituting Into (71c) and rearranging terms, give 
^33 ^3 " ^3 ( 1 ’*'^) 


where a^^ and are also defined. Next, from (71b) 
^2 f^^2“^23^3^ 


Its rounded value is 


^2 "" a«« C^ 2 "^ 23 X 3 (l+e)] (l+2e) 


3o 9 (l+2e) 


= Xg + 


2eb2 - 4e 333 X 3 
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X 2 “ Xg - 


Scbp * ^23X3 

“22 


A 

"X2 - 


2s bg - 4e 923^3/ 


( 76 ) 


Substituting (76) and (73) Into (71b) and rearranging terms, give 
^22^2 ^ ^^3 * ^2^^^^^^ 

A A A 

®22 ®23 h 

AAA 

Where a2g» 323* and b2 are also defined. Next, from (71a) 

xi =s {7 [^-*12X2-^13X3] 

Its rounded value Is 

X-j - — [b^(l+3e) - ^12^2^^^^^^ ” ai3X3(l+3e)] 

With the help of (72) and (75), 


(77) 


^1 ■*■ hIt r^^l *■ ^"^^ 12^2 *“ all ^2 **■ 


*12*23 


X3 - *13X3 


'11 '• ' ‘ “22 “ *22 
A A ^ 

Using the approximations X2 X2 and x^ ^ X3, x<j Is expressed in terms 

A A A 

of x^ , Xg, and X3 as 

a 

? i 

c^:^ c.{. 

(78) 


} 


J - js .., - A . * 


12 ..0 e . 4e /“12I23 . a ^0 \ 

‘22 13 ^^3/ 


Substituting (78), (76), and (73) Into (71a) and rearranging terms, give 
a^^ x^ + ai2(l+4e)x2 + ai3(l+3e) X3 == b-j(l+3e) (79) 


'n 


'12 


'13 
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A A A 


where ®i 3 » and are also defined. Finally, put (74), (77) 

and (79) Into a single matrix equation. 


nr^ll ®l2=ai2(l+4e) 'ai3*a,3(l+3eT] pH p^-b,(l+3c] 

0 ^22°®22 flgj-aggd+Se) Xg » 


822=822 


®33'®33 


A = A + M 


^ ^ + Ab 
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Its equivalent perturbation equation for evaluating rounding errors can 
be obtained by generalizing the result of (80), which 1s 



A Ik 


where 


“kj 

k«d 


®kj “ ^ 3yn'*'(n-j+3)E] 

l£k<j<n 

(85) 


otherwise 


b|^ = b|^[l+{n-k+l)E] 

k»l to n 

(86) 

Thus, 



[a|^.(n-j+3)e 

l£k<jin 


^""kj ° { 


(87) 

0 

otherwise 


% 

Ab|^ « b|^(n-k+l )e 

k « 1 to n 

(88) 
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11.4. A Numerical Example 

Consider Inverting the following matrix 



1.235 

-1.234 

3.25? 

A - 

1,023 

in 

CO 

CM 

in 

1 

0.921 


1.336 

2.120 

-8.235 


using computer of different finite decimal wordlength. Then evaluate 
the corresponding error norms using the procedure of Figure 1, The 
effectiveness of the procedure is examined by comparing these error 
norms to the corresponding actual error norms. The actual norms are 
approximately obtained by using a computer having a much longer decimal 
wordlength. The result is given in Table 1, which shows that error norms 
obtained by using the proposed procedure are indeed very conservative. 
Notice that error norm decreases with increasing wordlength. It is in- 
teresting to note that when the proposed method is used all error norms 
have the same mantissa. 


Table I 

Matrix Inversion Error Norms 


Wordlength: No. of places 
after decimal point 

Error norm 

By proposed method 

Actual value 

3 

3.706 X 10"^ 

7.782 X 10'® 

5 

31706 X 10"^ 

4.530 X 10"^ 

8 

3.706 X 10"^ 

4.672 X 10"^° 
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11,3. Resultant Rounding Error 1n the Inverse Matrix 

Two set of equivalent perturbations for A and ^ have been obtained 
to account for rounding errors. One set, A and ^ as given by (59) and (60), 
account for errors from triangularization. The second set, aA and A^as 
given by (89) and (90), account for errors from back substitution. The 
resultant equivalent perturbations for A and are given by the sums 

5A=’?J + aA (91) 


6^ = £ + Ajb ( 92 ) 

The resultant rounding error in x in the solution of 

A X = u,. j = 1 to n , (93) 

where u^j is the jth column vector of the 'tdentify matrix I, i-s given by 

Xj] (94) 

This equation is obtained in a way similar to that of (68). Rounding error 
*'-1 

in A is then given by 



The error norm of the computer's inverse matrix of A is therefore 

ll«(A-’)ll < ||A-’||.1|D|| (96) 

where the matrix 0 has been defined in (95). The relative error norm is 



(97) 


It is obvious that the evaluation of (96) or (97) involves a good 


deal of computation and should be done by a computer. Figure III-l Is a 


computation block diagram for this purpose. 
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The error bound 
for A"^ 


Figure Flow chart for evaluating rounding error 

bound of matrix Inverse. 
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12. Rounding Error Bound for Kalman Filtering 

Consider a process modeled by the following set of sequations. 
-k '•’k-1 ^-1 ° “ 

Z|( = H|^ X|^ + ^ dim ^ = tn 

E x(0) = Xg E[£(0) £(0)'''] = 
i!l,< 'V' N(0,Q|^) V|^ ■V. N(0,R|^) 

ECs!), v^] = 0 all j,k 

The Kalman Filter algorithm consists of the following equations. 


h. ^ *k-l ^-1 * ■ *^k 

^ (0) ~ ^ 

(98) 

Kk = Pkp ^kt^k^kp^k " Rk^’’ 

or = P^hJr-^1 

(99) 

^kp ° *k-l’’k-rk-l ^ '^k 


(100) 

^k “ ^kp “ *^p ^k ^kp 


(101) 


where x_ is the estimate of x. The present interact is to find the bound 
of the rounding error norm for ^ . For the sake of convenience, the 

A 

astrick will be dropped, and x.will denote the rounded value of the 
estimate. 
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Error in Rounded ,5^ 

’’kp '^k-l *’k-l '*’k-l ■*■ ^k-1 
Its rounded result 1s 

Kf * Qk-i 

= P|^ (l+5n+Te)-2nEq|^_, (102) 

The roundinq error is 

P|^p = (2n+l)eP|^p - 2neQ|^.^ i(2n+l)oPkp (103) 

Error in Rounded P,^ _ 

Pk = n - K^H^jP^p 
Its rounded value 1s 

P|^ = [I(l+o) - K|^H|^(l+m+tE)]P|^p(l+ne) 

= [I(l+n+le) - K|^H|^(l+n+iii+lE)]P|^ 

= [I(l+n+fE) - K|^H|^{l+n+iii+rE)][P|^p +. P,^p] 

^ P |5 + [{n+l)£l - K|^H|^(n+m+l)t]P|^p + Cl"K|^Hk^^kp 
= P|^[l+(n+m+l)E] - niE P,^p + [I-WP|,p ' 004) 

(h+pH-DeP^ -niEP^p^- [I-K,^H^]P^p , 

^ (n+m+l)ePj^ - mePj^p + [I-Kj^Pj^](2n+l ) Pj^p 

^||■ m ...^ l.ll■ ■ lll / 

(2n+l )eP|^ 

^ (3n+m+2)eP|^ • (105) 
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Error 1 n Rounded (99)^ that U 


Its rounded form Is 

1<l, = ■^I^Hj^n+n+nit;) 

By (104) and (103) 

P|^ (]+FwTe)Pj^ - RlePj^p + [I-Kj^Hj^](2n+l )eP|^p 

(2n+l)eP,^ 


" (l+3n+in+2e)Pj^ - niepj,p 
Assume Rj^ diagonal 

Than (106) becomes 

\ 1 C(l+3n+m+le)Pj^ - meP|^p]HjR|^^(l+e)(l+rFme) 
. = (l+4n+2m+3E)P^H[R'’ - mePi^pH^R"’ 


The rounding error is 

i^k i ('t"+2m+3)eK|^ - mePkpHjRt’ 


^ (4n+2m+3)eK|^ 


Error In Rounded (gg,^ is 

4 ^ ‘*'k-l’*k-l ■*■ ■ Vk-1^-1^ 

= [I-KpHp]iip.,Xp., + KpZp 


(106) 


(107) 


(108) 


(98) 
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F|^ = 

“ '^k ^ 

Then (98) becomes 

^ ° ’'k 4-1 !ik 

Develop the following rounded quantities. 

F)^ = [I(He) - K|^Hj^(l+ni^e)]<|>,^.^(l^ 

« [I(l+n7TE) - Kj^ ( 1 +4n+2m‘+3 e ) H,^ (1 +irHOTe ) ] <ji 

== Cl(l+n^e) “ Kj^H ( 1 +5n+3m+4e ) ] ^ j^_-| 

» F}.(H5n+3m+4e) - (4n+3n+3)el 

= K|^ ( 1 +4n+2m+3e ( 1 +me ) 

= K|^^(l+4n+3m+3e) 

= u^(l+4n+3m+3e) 

^ (l+5n+3m+4e)F|^ ><^...1 U|^(l+4n+3m+3e)(l+e) 

= (6n+3m+5e+l)F|^ + (4n+3i+4e+l )i^ 
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Define 

A 

F|^ = (l+6n+3in+5e)Fj^ 

^ = (l+4n+3ni+4e)Uj^ 

Then (114) can be written as 

4-1 " 4 

Case of k “ 3 Eq, (111) gives the exact as 

h "" ^ ^3^2-1 ^3^ ^ 

The rounded is 

A ^ ^ A /W ^ 

“3 ~ F2F2F^^(l+3n+3e) + F|jF2U;| (l+2n+3e) 


+ FgUgll+n+^e) + M3(He) 

Using (115) and (116), and combining terms, 

^ = F3F2F^^[l+3(7n+3m+6)e] + F3F2U^CHT8n+9m+17e] 

+ F3U^[l+nn+6m+l 1e] + i^3[1+4n+3m+5e] 

“ “ (3m+l )ep3F2Jii 

“ (10nH*3m+7)eF3]^ •- (17n+6in+l 3)ei^ 

The rounding error is 

= 3(7n+3m+6)e>^3 - (3n+l )eF3F2 U.-| 


015 ) 

(116) 


(117) 


(118) 


(1T9) 


- (1 On+3m+7)eF3 ^ - (17nt6m+13)eu^ 


( 120 ) 
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Its norm is bounded by 

I li 3 ll< 3 ( 7 n+ 3 m+ 6 )e|lx 3 l| - (3n+l )e| | F 3 F 2 M. 1 1 1 

- n0n+3m+7)e||F3U2l| - (17n+6m+13)E| IU 3 I | ( 121 ) 

Assume I |Fj^| | ^ F, | |U|^| \ ± U, and F^u, £ B for 1 » 0 to 2 , then 
11^311 £3 (7n+3m+6)e| IX 3 I I + (30n+9m+21 )eB 

1 3 | 3 ( 7 n+ 3 m+ 6 )||x 3 ||+[ 3 ( 3 n+l) + (7n+3nH-6)]Bl (122) 

Case^of k=4 Again, by (111), 

24 = ¥3>^2'^1 2o ^ ¥3'^2 iti + ¥3 Mg ^ >^4 24 023) 

Following similar derivation, gives the norm of rounding error as 

1 1^11 < B|4(7n+3m+b)|ix4|| + [4(3n+i) + (7n+3m+b)]B} (124) 

The general case_k From the equation pattern of (122) arid (124) for 
k=3 and 4, the general case is found to be 



(125) 


Eq. (125) is the main result of this chapter, which gives the bound 
of error norm for the rounded state estimate. The bound depends on the 
unit rounding error; k, the number of iteractif/ns; n and m, the dimension 
parameters of the process; |lX|^l|, the norm of the estimated state; and B, 
a quantity depends on Kj^, Hj^, and i^. The usefulness of this equation 
is at providing a general idea on the desired number of digits for the 
mantissa of the computer's floating number system. The following example 
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Exampl e . Consider a one year GP-B operation where relativistic data 
are taken every 10 seconds. Assume that the Kalman filtering involved in 
data reduction is also operated at 10 second iteractlon period. Then, at 
the end of one year period, the value of k would be k»365 x 24 x 3600/10 « 
3.1536 X 10 . Assume be a 10-vector and ^ be a 2-vector, so n»10 and 
m®2. Then, 

k(7n+3m+b) = 2.6490 x 10® 

k(3n+l) + 0.5k(k-l)(7n+3in+6) = 1.2949 x 10^® 

Just for the sake of discussion, assume only one term at the right-hand 
side of (125) dominates the result. If the first term dominates, one may 
estimate the desired 3 from 
k(7n+3m+6) 3=1 
so 

s = 0.3775 X 10‘® 

comparing to (8), the formula 3 = 0.5 x 10 , gives t _^9, therefore 9 

digits are desired for the mantissa of the floating point number system. 

On the other hand, if the second term of (125) dominates, one may estimate 
the desired 3 from 

[k(3n+D f 0.5k(k-l){7n+3m+6)] 3=1 
so 

6 = .7723 X 10"^® 

comparing to (8), gives t JL!17. Hence 17 digits are desired for the 
mantissa. 


G- 
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13. Remarks 


1. The main result of this chapter, given by (125), is a "first- 
cut" result which can probably be further refined to tighten the predicted 
bound while maintaining its reliability. This result was obtained after 
several different approaches to the problem had been attempted, 

2. It is desirable to find out the effectiveness of (125) ^y a com- 
puter emulation of the GP-B data reduction. This has not been done for three 
reasons. First, there is no time and resource allocated in the contract 
period covered by this report. Secondly, this emulation work is a major 
activity not foreseen at the start of the contract, therefore it was not 
planned. Finally, the computation details of the GP-B data reduction, 

being developed by the GP-B group at the Stanford University, is not 
completely available to the present contract. 

3. Besides Kalman filtering the GP-B data reduction also contains 
other types of computations which should also be included in the rounding 
error analysis, Attention to a more complete error analysis is intended 
for for the next contract period. 

4. During the present contract period only the finite wordlength 
Induced computation errors were investigated. The truncation error in- 
volved in the chosen method of numerical integration for the GP-B data 
reduction has not been studied, thus its improtance in the overall error 
analysis is not known at this time. The effect of Truncation error will 
be investigated in the future. 


5. Other indirect errors exists, which also enter the result through 
GP-B data reduction. These errors are caused by Imperfection in compensa- 
tion data. Effect of this type of errors has not been Invtistigated, but Is 
Intended for future study. 
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CHAPTER ly 


MEASUREMENT GEOMETRY FOR GP-G EXPERIMENT 

This chapter describes in detail the development of a measurement 
model for the GP-B experiment. The model represents the relationships 
among the orientations of various components of the experiment, including 
the directions of the reference star, gyro spin axis, readout normal, and 
the spacecraft’s roll axis. With this model, the effect of spacecraft 
rolling and of component misalignments can be investigated through analysis 
and computer simulation. 

For the convenience of analysis, six coordinate frames are defined 
with the help of Figure IV-1. All frames are co-origined at the spacecraft 
center of rotation. 

1. Absolute frame [a-frame: x^, y^, — This frame is 

stationary in orientation with respect to the universe, 
i.e. With respect to the ideal distant stars. Its Xg^-axis 
is along the line-of-sight (LOS) to the reference star 
Rigel. 

2. Intermediate frame [i -frame: x^, y^, — This frame is 

stationary in orientation with respect to the universe, 

Its x-| axis is along the roll axis of the spacecraft. This 
frame is related to the absolute frame by Euler angles 
«g, and oij* 

3. Roll frame [(fi-frame: x^, y^, z^] — - This frame is fixed to 
the spacecraft with its x^ axis along the roll axis of the 
spacecraft. Therefore x^j^ and x^ axes coincide. 
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Figure IV-1. Definition of directions. 


4 . Telescope frame [t-frame: x^, y^, z ^]‘ — This frame is 
also fixed to the spacecraft, with its x^ axis along the 
boresight of the telescope. The misalignment between this 
frame and the roll frame is represented by Euler angles 

3 ] 9 ^2 » and • 

5 . Readout frame [r-frame: Xy,, y^,, z^,] — This frame is again 

fixed to the spacecraft, which represents the orientation 

of the gyro readout ring* Its axis is normal to the 
plane of readout ring. The misalignment between this frame 
and the telescope frame is represented by Euler angles 
Y’]> Y2* ^3* 

6. Gyro frame [g-frame: Xg, yg, Zg] — THis frame is 
inertiaVly fixed, but not absolutely fixed. The relative 
angular motion of this frame with respect to the absolute 
frame is the relativistic drift of the inertial space 
around the GP-B gyro. The orientation of this frame with 
respect to the readout frame is represented by Euler angles 
6-j, 62, and iSj# while that with respect to the intermediate 
frame is represented by Euler angles e-j, Bg, and 63. 

Relationships among all frames are also depicted in Figure IV- 2 . Ideally, 
3's, y's, and 6's are zero, meaning that the roll, telescope, and readout 
frames are all lined up. The nonzero values for a's represent the 
deviation of telescope's boresight from the LOS to' star. The e's will 
not be zero, they are the relativistic drift angles. It is assumed that 
angles, a's, 3's, y's, 6's, and e's are all sufficiently small that 
small angle approximations of trigonometry apply. 
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Figure IV-2. Relationships among coordinate frames. 
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In the experiment, an optical system meaures the angle between the 
LOS to star and the telescope's boresight. This amounts to measuring 
the projection of x^^ axis of the absolute frame on the telescope frame, 
With all misalignments included, this relationship 1s modeled by the 
following matrix equation. 
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The gyro and its readout ring measure the relative angular displacement 
between the gyro spin axis Xg and the readout normal z^. This amounts to 
measureing the projection of Xg on the readout frame. Note that the readout 
frame rolls with the spacecraft since it is body fixed. With all mis- 
alignment included, the projection of Xg on (x^., y^, z^) Is modeled by 
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Gyro spin 
axis unit 
vector in 
r-frame 




Transformation to 
r-frame. Accounts 
for misalignments 

Y-j » Yg* Y3‘ 


Transformation to 
frame. Accounts 
for rolling about 
X- axis 


Gyro spin axis 
unit vector 
‘in g -frame 

( 2 ) 


Under the ideal condition, 3's and y's are all zero. Then (1) and (2) 
reduct to 




Solving (3) and (4) for a‘s and e's, gives 



offset of star's LOS 
from x^ axis 


(4) 


( 5 ) 
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( 6 ) 


Notice that when small angle approximations are used a-j and e-j cannot be 
determined. Since x^ and x^, axes are col inear, subtracting (5) from (6) 
yields the gyro's relativistic drift angles iji's with respect to the inter- 
mediate frame. Because of the absence of , the transformation from the 
intermediate frame to the absolute frame is simply via an identity matrix. 
Therefore the relativistic drift anglesof the gyro with respect to the 
obsolute frame are given by the same ij^'s. 




! — — 


’I's 


; “3"®3 


ro 

L 


_®2-“2 



Cwt -Sojt 
Swt Cwt 






( 7 ) 
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The Q ' matrix in (7) represents the deroll operation to bb done by the 
data reduction computer. 

To study the effect of misalignments on the accuracy of the experi- 
ments, (1) and (2) are used. In there, 3's and y‘s are unknown quantities, 
since the effect of known misalignment can be compensated computationally. 
Angles a's and e's are desired, and they are estimated from the measurement 
by the deroll computation, giving 
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The errors due to misalignments are 
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The effect on the relativistic drift angle determination is 


( 10 ) 


( 11 ) 
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This completes the present development of the measurement model for 
the GP-B experiment. The model will be used for investigating the quanti- 
tative effect of various misalignments on the experiment data. 
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